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The boundary of the first order chiral phase transition region is studied in the m-^ — rriK — /is 
space using the one-loop optimized perturbation theory for the resummation of the perturbative 
series. Chiral perturbation theory for mesons and baryons is used for the T = parametrization 
of the model. The surface of second order transition points bends with increasing /is towards the 
physical point of the m-„ — mx mass plane allowing for the existence of the critical end point in the 
fJ-B ~ T plane at the physical point. The location and scaling region of the CEP is explored. 
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In the past few years much effort was invested in the mapping of the phase diagram of strongly interacting matter. 
Apart from the theoretical importance of the phase transition for understanding of the physics of the early universe 
, the continuous interest in this issue is also triggered by the fact that some parts of the phase diagram are within the 
—i ' reach of current and future heavy ion collision experiments. In QCD and also in its effective models the parameter 
K*" , space contains both experimentally tunable quantities like the temperature and various chemical potentials (isospin, 
2? ' strange, baryonic) and parameters like the quark masses which have a fixed value in nature (physical point). At zero 
chemical potential the nature of the phase transition at the physical point was established only very recently. The 
continuum extrapolation in the lattice investigation of Refs. [l|, 0] using staggered fermions showed that the transition 
is of analytic cross-over type. How far the physical point is located from the boundary of the first order transition 
• region, expected on theoretical grounds in the limit of vanishing quark masses, is not yet definitively quantified. Both 
lattice simulations with improved staggered quarks (p4 action) ]3| and effective model studies [1, [^, [a, 0| show that for 
degenerate quark masses the first order transition region extends to pseudo-Goldstone meson masses with values well 
' rS i' below the physical mass of the pion. 
I [ For non-vanishing baryonic chemical potential the first order transition region persists and its second order boundary 
forms a critical surface in the niuM — rus — space on which the transition as a function of temperature is of second 
(— I \ order. Generally it is expected that for high values of the chemical potential the phase transition is of first order, so that 
• • . the critical surface bends towards the physical point (positive curvature) . In this case at physical quark masses there 
^ ' is a critical end point (CEP) in the chemical potential-temperature plane at which the first order transition line ends. 
In this plane, for chemical potential values smaller than the value at CEP the phase transition changes to a cross-over. 
The CEP was found in staggered lattice QCD first for a non-physical value of rriud [1| and then for physical values of the 
quark masses [§| (pion to rho meson masses tuned to its physical value) . The location of the CEP was also estimated 
in Ref. using Taylor expansion around /is = 0. However, the question of the existence and location of the CEP is 
not yet settled for two reasons. First the continuum extrapolation is not yet performed and not all the lattice methods 
used to address the notoriously difficult problem of performing simulations at finite chemical potential could find the 
CEP. According to the simulation of Ref. [ll| performed with imaginary chemical potential the critical surface turns 
out to be extremely quark mass sensitive and contrary to the standard expectation seems to have negative curvature 
(the first order transition region shrinks with increasing ^s). This would imply the absence of the CEP in the vaHdity 
range of the method used in Ref. estimated to be < 500 MeV, but leave the possibility of its appearance above 
this limiting value. 

In this paper our aim is to study the phase boundary of the first order transition region of the chiral phase transition 
in the TOtt — ttik — f^B space with special interest for the location of the critical end point in the /is — T-plane occurring 
at the physical point of the mass-plane. For this investigation the SU{3)l x SU{3)r chiral quark model is used as 
an effective model of strong interactions. This model has the same global symmetry features as the QCD. In view of 
the result of which shades doubts on the existence of the CEP, we investigate thoroughly, by studying the predicted 
mass spectrum, how reliable is the existence of the CEP predicted by our model (Sec. [ill- Then by continuing the 
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parameters of the model as functions of 171^^ and rriK using the chiral perturbation theory (CHPT) for meson and baryons 
we study the shape of the critical surface which bounds the region of chiral first order transitions in the rriTr — rriK — ^J.B 
space (Sec. lIII A|) . For physical values of TOtt and ttik the location of the CEP in the ^jlb — T plane will be determined 
in Sec. HiTbI We conclude in Sec. JVl 



II. THE MODEL AND THE DETERMINATION OF ITS PARAMETERS 

The Lagrangian of the SU{3)l x SU{3)ii symmetric chiral quark model containing explicit symmetry breaking terms 
is given by 

L ^ ^Tr {df^M^d'^M + ttiImHI) - /i (Tr [M'^ M)f - /aTr {M^lf - g (det(M) + dct(Mt)) + e^a^ + egag 

+^(^^-(7i.Af5)V'■ (1) 

The quark field '4) = (w, d, s) has implicit Dirac and color indices. The two 3x3 complex matrices are defined in 
terms of the scalar and pseudoscalar tt^ fields as M = Y^^i=o{'^i + 1^^%)^ and A/5 = X]^=o l^'^i + ^l5'^i)^ij with 

Xi : i = 1 . . .8 the Gell-Mann matrices and Aq := y^l- From the original 0-8 basis we switch to the non-strange (x) 
- strange (y) basis by performing the orthogonal transformation 



V3 V 1 -V2y ^"8 



(2) 

applied to the scalar, pseudoscalar mesons and external fields, that is u e {cr, tt, e}. Then one has 



1 ^ 1 
M5 = — = V'lcr, + ij5T^i)Xi + — =diag(cra: + il5'!^x,o-x + ij5T^x,V2{ay + 175%)), (3) 
V 2 ,=1 V2 

and similarly for M. 

In the broken phase the fields are shifted with their expectation values x := (ax) and y := {<Jy), which are obtained 
from (tTo) and (cts) cf. Eq. No isospin breaking is considered which means that the u and d constituent quarks are 
degenerate in mass: = Md = gFx/2. The mass of the strange constituent quark is = gpy/^f^- The explicit 
expression for the mesonic sector of the Lagrangian obtained after performing the traces can be found in [3, The 
sum of terms linear in the fiuctuations gives two equations of state which at one-loop level are 

= (— ) = -ex - mix + 2gxy + if^xy^ + 2(2/i + /a)^^ + ^ ^ {a,a,) + ^{{uu) + {dd)), (4) 

= = -cy - mly + gx' + if^x^y + 4(/i + f2)y' + J^ tl^ ^ (a^a,) + ^{ss), (5) 

where a goes over the pseudoscalar and scalar fields and z e {1 . . .7,x,y}. The non-zero three-point couplings t%. . and 
t^. . are given in Appendix C of Ref. 0|. Note, that in a given multiplet the three-point couplings are degenerate, e.g. 
til — ^22 = ^33 • The coefficients of the quadratic terms in the fiuctuations are the tree- level masses (see Table 1 of 0). 
The mixing in the x — y (0—8) sector is refiected by the fact that the off-diagonal xy component of the mass matrix 
is non- vanishing. In Appendix [A] we sketch how the propagators of the mass eigenvalues (77,77' and a, /o) are obtained. 

In Ref. 01 a method was proposed to determine the parameters of the Lagrangian. At the core of this procedure there 
is a resummation performed using the so-called Optimized Perturbation Theory (OPT) of Ref. in order to avoid the 
appearance of negative propagator mass squares in the finite temperature calculations with one-loop accuracy. In the 
OPT the mass parameter — ttiq of the Lagrangian, which could be negative if the model exhibits symmetry breaking, 
is replaced by an effective (temperature-dependent) mass parameter 7n^: 

imass = ^m'^TrM^M - ^(771^ + m'^)TvM''M = ^m^TvA-I^'M ~ ^Am^TTlVpM. (6) 

771^ is determined using the criterion of fastest apparent convergence (FAC) and the finite counterterm Am'^ is taken 
into account first at one-loop level. 

Compared to apart from the fact that we have to fix the Yukawa coupling g^, we have changed the parametrization 
in two aspects. 1. In place of the one- loop relations of the Partially Conserved Axial- Vector Current (PC AC) we 
use only the tree-level relation for the pion decay constant, because it turned out that in contrast to the one loop- 
bosonic contributions, which are explicitly renormalization scale independent, the fermionic contribution depends on 
the fermionic renormalization scale. 2. We have modified the gap equation for the pions and defined the one-loop 
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pion mass as —iG^^{p = 0) = instead of defining it as a pole-mass. This is because previous investigations [tI. fl5| 
showed that an unavoidable feature of the OPT method is that at finite temperature the solutions of the gap equation 
and of the equations of state cease to exist above a certain temperature, unless we impose arbitrarily that the finite 
temperature part of the mesonic bubble integral is taken at vanishing external momentum. 

In what follows we present the equations used to determine at T = the 9 parameters of the Lagrangian: the couplings 
™0i /ii /21 9i 9f, the condensates x, y and the external fields Cy The physical content of the one-loop pseudoscalar 
self-energies are shown diagrammatically in Fig. 1 of here we give expHcitly only the fermionic contributions. The 
integrals are given in Appendix [b1 The renormalization of the mesonic part was discussed in where the counterterms 
i5mQ, 5g, Sfi and 6/2 were given by Eqs. (4)-(7). Since the fermions are taken into account perturbatively one can 
easily see from Eqs. Q and ^ that using cut-off regularization the renormalization of the momentum- independent 
fermionic parts can be achieved with the following counterterms 
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A is the 3d cut-off. The renormaHzation of the momentum-dependent fermionic parts can be achieved with a wave 
function renormalization constant in the mesonic sector: 5Z = — tIt- In 4-- 

The FAC criterion used to determine the effective mass is implemented by requiring that the one-loop pion mass 

= -ml + (4/1 + 2f2)x^ + 4/iy2 + 2gy + = 0, m,, M„) (8) 

I 

stays equal to the tree- level pion mass (M^ ^ TOtt)- The fermionic contribution in the self-energy = 0) is 

= 0) = -2g2 T;^(M„), (9) 

where M„ = gpx/2, and Tp{Mu) the fermion tadpole defined in Appendix [Bl 

Using the expression of the tree-level pion mass the equation above results in a "gap" equation for the effective mass 

= -ml + = 0, m,(m2), M„). (10) 

We introduce the expression above for the effective mass in the expression of the tree- level pion mass. For all the other 
tree-level propagator masses appearing in the self-energy we replace the effective mass with the pion's mass using 
its tree-level expression. Thus we arrive at the following gap-equation for the pion mass 

ml = -ml + (4/1 + 2f2)x^ + 4/iy2 + 2gy + S^b = 0, m,(m,), M„). (11) 

This is the first equation from the set of equations used for the parametrization at T = and has a distinctive role in 
our investigation because in the thermodynamical calculation one has to solve its T-dependent counterpart in order to 
determine m^(T). 

The next 2 equations require that the one-loop pole mass of the kaon and eta are equal to the corresponding physical 
masses. The one- loop equation for the kaon and eta pole masses read (see Q for their derivation): 

Ml = -ml + 2(2/1 + f2)ix^ + y') + 2/2z/ - V2x{2f2y - g) + Re {^^(p^ = M|., m,)} , (12) 
^Re |m^^^ + S^^^ (p2 = M^,mi) + m^^^ + S^^^ (p^ = M^, m,) 
- ^(m2^^ + E,^^ (p2 ^ Ml m,) - m^^^ - {p^ = M^ m,))^ + 4{ml^ + E,^^ {p^ = M^ m^w] . (13) 

The fermionic contribution to the self-energies of the kaon and of the different components in the x — y mixing sector 
are: 

^xip) = -9I [Tf(M«) + Tp{Ms) - {p" - (Mu - Msf)lF{p, M„, Ms)] , (14) 
S^.. {P) = -9I [2Tf{Mu) - p'If{p, MJ] , S^^^ (p) = -g% [2Tf{Ms) - p'lpip, M^)] , [p) = 0, (15) 

where the fermionic bubble integral If is given in Appendix [b1 We implemented the FAC principle for the kaon by 
requiring additionally S(p^ = ^ix) — 0. The gap equation for the pion together with these three equations are sufficient 
to determine mg, /i, /2 and g, if we know x, y and gp- The tree-level PC AC equation fixes x = Jt^. Then the Yukawa 
coupling is obtained from the tree-level relation gp = 2Mu/x, where M„ is the non-strange constituent quark mass. 
With the value of gp fixed we determine y from the strange constituent quark mass y = y/2Ms/ gp- The values of the 
external fields are determined from the two equations of state ^ and ([5]). 

We are forced to use the tree-level quark masses because the one-loop quantum correction to it turned out to be very 
large. This could be kept reasonably small only by using unrealistically low mass values for the tree-level constituent 
quarks. 
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FIG. 1: Contour levels for the quantity R (Eq. iflG]! ^ which gives the average percentage for the deviation of the predicted 
spectrum from the physical one. Shown is also the order of the phase transition on the T = 0//iB =0 axes of the T — /is-plane 
as a function of the fermionic and bosonic renormalization scale // and It (c stands for cross-over). 



Since we are working with one-loop self-energies, two renormalization scales emerge, one for the renormalization 
of bosonic integrals, h and one for fermionic integrals. If. The parameters were determined as functions of these 
renormalization scales h and If using as input the values of the physical quantities = 138 MeV, niK = 495.6 MeV, 
rrir, = 547.8 MeV, = 93 MeV, A/„ = 313 MeV and Ms = 530 MeV. The values of the constituent quark masses are 
related to the masses of the baryon octet components as M„ « and « (Ma + Ms)/2 — 2M„. Constraints 

on the values of h and I / can be obtained by confronting the predicted part of the tree-level and one- loop level mass 
spectrum with the physical data. To quantify the deviation a quantity is introduced such to reflect the best the spirit 
of our parametrization 

1 Kee_^phy.| ^ ,^tree _ ^1-loop | 

\T\^ mfy^ ' 

where T = {77, 77', ao, /o, cr}, L = {?/, aq, k, /q} and \T\ ~ 5, \L\ ~ 4. Since the propagator masses that enter the one-loop 
integrals and determine the quantum corrections are the tree-level masses we measure their average deviation from a 
given physical spectrum. We also included into the test quantity, R the average deviation of the one-loop masses from 
the tree-level ones. The one- loop masses are determined from the corresponding propagators through the equation 
'ReiG~^{p^ = "m.^ i_\oop) = 0. There are two notable omissions from the right hand side of Eq. (fT6| . to^"'°°p is omitted 
because in a relatively large range of the If — /^-space to be scanned ReiG~^(p^ = Jti^ i-ioop) ~ solution. 
m^i^'^'^ is omitted because it is actually independent of lb and If. This becomes evident if we use the tree- level relation 
tree = ('^i'/if ~ '"^^ f ir) / {f K — fir) where at tree-level fK = x/2 + y/\/2 and according to the parametrization 
above neither x nor y depend on 4 and If ^. The pseudoscalar physical spectrum used as a reference is completed 
by TO,,' = 958 MeV. As for the scalars, it is still not clear which mesons form the lightest nonet. Possible candidates 
are the isoscalar a and /o(980), the isovector ao(980) and the isospinor k(900). For the isoscalar masses we choose 
TOct = 700 MeV and to/q = 1370 MeV, because this parametrization of the model cannot accommodate an /o lighter 
than 1 GeV. 

The contour plot of the average percentage for the deviation from the physical spectrum measured by R can be seen 
in Fig. [TJ Based on this plot we have chosen the point marked with a star in the figure for which h = 520 MeV and 
If — 1210 MeV. At this point the 1-loop sigma mass can be obtained using the pole-mass definition given above and 
one obtains mo- = 614.2 MeV. The one-loop mass of /o is to/„ — 1210.9 MeV. Around this optimal point the quantity 
R can be completed by the contribution of the one-loop mass of the a meson and the new analysis confirms the choice 
of Ib and If given above. 

As indicated in Fig. [T] the order of the phase transition on the axes of the — T plane turns out to be sensitive to 
the value of the renormalization scales. It is very interesting to note, that the requirement of minimal deviation from 
the physical spectrum results in values of the renormalization scales for which the transition is of first order on the 
T = axis and of cross-over nature on the = axis of the fiB — T phase diagram. This feature turned out to be 
independent of the way we quantify this deviation, that is the way we define the quantity R. In conclusion, the most 
reliable parametrization of this model positively predicts the existence of the CEP. 



^ Note, that with this parametrization the kaon decay constant turns out to be /k = (1 + Ms/Mu)fn/2 = 125.23 MeV which is 10% larger 
than the physical value 
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III. THE PHASE DIAGRAM IN THE m„ d - m» ^ib SPACE 

In order to determine the order of the phase transition in the m„ ^ — rus — hb space we have to determine the 
temperature dependence of the order parameters. This is obtained by solving the two equations of state and the gap 
equation which form a set of coupled equations: 

= -e, + m^x + 2gxy + ifixy^ + (4/i + 2f2)x^ + ^ J,t^T^{m,{m^)) - igpMuT^iMu) + Am^x , (17) 

i 

= -ey + mhj + gx^ + Af^x^y + 4(/i + /a)^^ + ^ J^t^T^{m,{m^)) - 2V2.gFM,T^(M,) + Am^y , (18) 

i 

ml = -m2 + (4/i + 2/2):r2+4/iy2 + 2gy + E,(p = 0,m,(TO^),M„). (19) 

Ji stand for the isospin multiplicity factors which are Jtt.oo = 3, Jk,k — 4, and J-q.^y .ajo — 1- The explicit renormalized 
expressions of the temperature dependent integrals pijn) are given in Appendix [BJ In the notation above we made 
explicit that the tree-level bosonic masses depend on the pion mass determined by the gap-equation. A first order phase 
transition is recognized by the multivaluedness in a given temperature and/or chemical potential range of either of the 
two vacuum expectation values. The point where a first order phase transition goes over into a cross-over is identified 
as a second order phase transition point. 



A. The surface of second order phase transition points 

In order to obtain the surface of the second order phase transition points in the m,j ^ — — hb space we have to 
determine the variation of the parameters of the Lagrangian with the pion and 1-caon masses. In addition to the mass 
of eta and the decay constant of the pion, the constituent quark masses are also used in the parametrization, so their 
and rriK dependence is required as well. For and m,, we use the formulas of the CHPT in the large Nc limit at 
0{l/p) as in Ref. 0| 

U = / (^1 + 4^5^), (20) 

^1 = ^ + y(2L8-L5) ^ ^' ■ (21) 

The values of the constants are: L5 = 2.0152 • 10"^, Lg = 8.472 • 10"^ and / = 91.32 MeV. These quantities are 
independent of tti^ and mx, and were determined in Ref. Q from the condition that at the physical point the pion 
and kaon decay constants and 777^ take their physical values. 

In a simple constituent quark picture the constituent quark masses are related to the values of the baryon masses 
whose dependence on 777^ and rriK is obtained using the formulas of the CHPT for baryons at 0{q^) [q denotes 
the momentum) . To this order in the chiral expansion the masses of the baryon octet are given by 

Mb = Mo - 2feo(m2 2 + 7772^^2) + bole + &f7b - + a%m], + alm^^ , (22) 

where B e {iV, S, A, S, }. a^, a^, are simple expressions of two low energy constants D and F, the coefficients 
characterizing the lowest order baryon Lagrangian (see Eq. (6.9a) of Ref. for the expHcit expressions). The 
expressions of 7^ and 7^^ read as: 

7^ = 7h = -4m2, 2, 7^ = -7l = 4(777^,2 - <2), lE = -^m^^^, ^ -^rnl, 7^=7^ = 0. (23) 

The subscript 2 refers to the fact that the squared meson mass is taken at the lowest, 0{q^) order. Originally these 
lowest-order mass expressions appear also in the 0{f~'^) term, but it is allowed to replace them with their physical 
values. For the O(f^) piece we use 777!^^ 2 = "^1/ ~ ^^^^ /-2(2L8 - L5) for M = TT,K and Eq. ^ to relate the 
lowest-order mass with the physical one. Note, that due to the lowest order Gell-Mann-Okubo relation the first term 
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in Eq. ^ is 

The unknown constants are Mq, the baryon mass in the chiral Hmit and 60, bo, bp the coefficients of the lowest 
order symmetry breaking part (contact contribution). Since Mq and 60 appear in the same combination for all the 
masses an additional relation is needed in order to disentangle them. This relation is usually chosen to be the so-called 
nucleon sigma-term defined in terms of the proton matrix element of the 11 component of the sigma commutator 
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The surface of second order phase transition points. 



{p\(Tii{0)\p) /2Mp, where Mp is the proton mass. ^ At 0{q'^) the sigma-term reads 



m: 
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2{2bo + bo + bF)ml 



647r/2 



'TT 



(24) 



The four constants are determined from the expressions of Mjv, Afc, (Ms + M\)/2 and cttttv using as input D = 0.8, 
F = 0.5, A/at = 938.92 MeV, Ms = 1193.15 MeV, Ma = 1115.683 MeV, M^ = 1317.915 MeV and a^^ = 45 MeV. The 
obtained values are: Mq = 939.82 MeV, 6o = -0.869 GeV-\ bo = 0.0363 GeV"! and bp = -0.582 GeV"^. We keep 
these values fixed when determining the constituent quark masses on the m^r — niK mass plane as Mu — Mn /i and 
Ms = (A/a + My,) 12 — 2M„. Then at fixed values of the renormalization scales the parametrization described in secUll 
is performed for each point of the mass plane. 

The result obtained after finding the nature of the phase transition by solving Eqs. (fT7|) . (fTSj) . ifTO ]) is displayed in 
Fig. [21 The second order surface is shown only in a restricted region of the mass plane for two reasons. On the one 
hand the parametrization breaks down close to the diagonal of the mass-plane. Along the diagonal the masses become 
degenerate and we cannot determine the same number of parameters from fewer equations. On the other hand, for 
large values of the kaon masses to_r- > 400 MeV we see that the surface bends away from the to^ = axis (apparently 
no crossing occurs). This is an unphysical behavior since within the linear sigma model, that is without fermions, it was 
shown in Ref. 01 that for /is = the boundary of the first order transition region approaches the = axis, allowing 
to locate there a tricritical point (TCP). This anomalous behavior originates from the way the baryon masses depend 
on rrLK'- all of them start to decrease with increasing rriK, for niN at around rriK = 300 MeV, for toa at rriK ~ 400 MeV 
and for ms and this occurs at around uik = 500 MeV. This means that above uik = 400 MeV one can not trust 
the 0{q^) formulas of the CHPT for the baryon masses. Based on the uik dependence of the baryon masses presented 
in we can say that had we used the more compHcated 0{q'^) formulas we would have observed the bending away 
of the surface from the = axis at around niK = 600 MeV. The result of Ref. Q shows that even at this value of 
the kaon mass we would be far from the scaling region of the TCP. 

We can see that the surface grows out perpendicularly from the mass plain {^b = 0). This means that the critical 
points iig{mTr,mK) which are close to the mass plane are extremely sensitive to the values of to^ and ttik- This was 
observed in lattice simulation using imaginary chemical potential, where this is the region which is traversed when 
performing the analytical continuation to real chemical potential We can also see that the tangent plane to the 
surface has a decreasing angle with the mass plane as one approaches the critical point which corresponds to the physical 
masses. Here the dependence of ^^(m^, m^) on the masses is milder. 



Now we turn back to the physical point in order to study the location of the CEP in the — T-plane. There are 
many effective model studies published in the literature (see Ref. [l^ for an extensive list), the major part of them was 
done with two fiavors. 



^ The sigma commutator is defined as a''''{x) = [QJ^(a;o), [Q^i^o) , Hshi^)]] , where Q^(xo) is the axial charge operator and Hsh = qMq is 
the chiral symmetry breaking term of the QCD Hamiltonian (Ai = diag(m„ , rrid, ms)). 



B. The CEP at the physical point of the mass-plane 
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FIG. 3: The phase diagram corresponding to the point denoted with a star in Fig.[T] Shown are the universal chemical freeze-out 
line from Ref. [l^l and the CEP obtained in the lattice The arrow at the CEP shows the temperature direction of the 3d 
Ising model. 

As we have discussed in Sec. [ill reasonable values, namely If = 520 MeV and lb = 1210 MeV have been chosen 
for the renormalization scales. In the cross-over region we determined for increasing values oi fJ-B the pseudo-critical 
temperature Tc through the location of the peak of the derivative of the non-strange order parameter x with respect 
to T. The peak diverges at the CEP. For higher values oi in the first order transition region the spinodals were 
defined as the two turning points of the two- valued order parameter x as function of /is at fixed T. Since we have not 
calculated the effective potential, we defined the first order line through the points between the two spinodals, where 
x{iib) has an inflection point at flxed T. This line together with the two spinodals, plus the cross-over transition line 
and the location of the CEP (Tcep = 74.83 MeV, ^b^cep = 895.38 MeV) can be seen in Fig.H The fact that the CEP 
value of the temperature is smaller and that of the //s is larger than what is found on the lattice (Tcep = 162 ± 2 MeV 
and /J.B,cEP = 360 ± 40 MeV, cf. seems to be a common feature of all results obtained in the linear sigma model 
and the Nambu-Jona-Lasinio model. 

In Fig. [3] we display also the universal chemical freeze-out curve of Ref. [131 which is parametrized as T{^b) = 
a-bfi%- cfi%, where a = 0.166 GeV, b = 0.139 GeV"^ and c = 0.053 GeV^ We can see that for > 400 MeV 
our phase transition line lies above the curve which separates a region dominated by inelastic processes (high-T region) 
from a region dominated by elastic processes. The parameters of the freeze-out curve {iib and T) reproduce in thermal 
models the particle yields measured in experiments performed at different beam energies. It is called universal because 
is obtained from the condition that the average energy per average number of hadrons is approximately 1 GeV. The 
sensitivity of this curve to other freeze-out conditions is studied in Ref. [l^l • It is expected that the chemical freeze- 
out takes place below the true phase transition line. In the resonant gas model with rescaled mass spectrum such 
as to reproduce TOtt = 770 MeV of the lattice simulation with 2 flavor QCD using Taylor expansion of the fermion 
determinant, the line of fixed energy density reproduces the transition line of the lattice. The set of physical resonances 
expected to appear in the (2+l)-flavor QCD results in a line of flxed energy density which is close to the chemical 
freeze-out line in the meson dominated region (/is < 400 MeV) and lies above it in the baryon dominated region 
{HB> 400 MeV) 

There are some quantities like the width of the peak of the chiral susceptibility and the curvature of the cross-over 
line at hb = 0, which can be compared with their values measured on the lattice. The width of the chiral susceptibility 
as a function of T gives information on the strength of the cross-over transition. The continuum result of Ref. [l| is 
ATc{x-4}ip) = 28(5)(1) MeV. We obtain ATc{xx) = 15.5 MeV which means that the transition is more abrupt in our 
case ^ . It is interesting to remark that for the simulation of Ref. [§| in which the CEP was found the width is about 
an order of magnitude smaller than in the continuum limit [l^ . Then a natural expectation is that in the continuum 
limit the location of the CEP would move to higher values of fiB, that is a larger value of ^b would be required to turn 
the phase transition into a first order type if one starts with a week cross-over at /is = 0. 

Another quantity which can be compared is the curvature at //s = which is measured by C := "^c^^I^^^q- Our 
value C = -0.09 is 35 - 40% larger that the lattice result C = -0.058(2) for Nf = 2 + 1 of Ref. l3ur pseudo- 



^ It is worth to note, that the susceptibility which can be defined in this model, namely x = dx/de^ can be connected with the chiral 
susceptibility of the light quarks X.^^ = d{'ipil)) / drriu ■ This is based on the following relations: (i/ii/') ~ = Borriu (lowest order 

relation in CHPT) and = tx/x (Goldstone theorem). From these immediately follows: ~ ^X- 
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FIG. 4: The chiral susceptibility as a function of the distance from the CEP. For the sake of visibility, the lines from left to 
right were shifted to the left by 0, 10° \ 10" ^ lO" ^' and lO" *, respectively. 

critical temperature at vanishing chemical potential Tc ~ 154.84 MeV is quite close to = 164(2) MeV obtained 
in Ref. [§| but is even closer to what was obtained in the continuum limit in Ref. namely Tc = 151(3) MeV. A 
collection of curvatures obtained in different lattice approaches and for different flavor numbers is given in Ref. (2ol |. 
where it was observed that the curvature is considerably decreased if the coefficient g of the determinant term in the 
Lagrangian decreases with the temperature. The rapid decrease of the anomaly near the chiral transition temperature 
is indicated by lattice simulations. In NJL model the T-dependence of g can be extracted by fittin g th e lattice result 
on the topological susceptibility with an analytic formula of the susceptibility . Following Ref. |22l | we considered 
g{T) = exp(— (T/Tq)^). As a consequence the curvature was reduced to C = —0.08. This means that the restoration 
of the Ua{^) symmetry has an influence on the shape of the crossover transition line and on the location of CEP. More 
detailed dependence of g on T and /xs is needed in order to assess this influence on the location of the CEP. 

We have studied the shape of the critical region based on the chiral susceptibility x = dx/dcx- This is important in 
phenomenology, since according to Ref. [2BI if the critical region is larger that the region where the so-called focusing 
effect (see Ref. [131) is realized, then the interactions cannot be neglected and the use of resonance gas model could 
be questionable. We found that the critical region is heavily stretched in the direction of crossover transitions line as 
shown in the smaller picture on FiglSj in which we depicted the ratio x/x^'^*'", where x^'"'"' is the chiral susceptibility 
of a free massless quark gas, given by x'^'"*''' = T^/G + (IStt^). The contours shown are 1.2, 0.8, 0.6 from inside out. 
The highest values are concentrated around the pseudo-critical line in a 1 — 2 MeV wide elongated region. The shape 
of the critical region is similar to what was observed based on quark number susceptibility in Ref. [131 in the two-flavor 
effective QCD. 

We determined the mapping of the temperature axis of the Ising model onto the ^b—T plane at the CEP which is in 
the universality class of a 3d Ising model. To do this, we have measured, through the relation x ~ (IM-B — M-B,cep| cosa-|- 
\T — Tq^p \ sina)"'*', the critical exponent of the chiral susceptibility along different directions pointing toward the CEP 
and characterized by an angle a € [0, 360°] measured from the positive axis. Going with increasing /is along the 
direction tangent to the cross-over line at the CEP having a = 320.85° we obtained 7 = 1.01, which corresponds to the 
mean-field exponent 1. This direction corresponds to the temperature axis of the Ising model and it is shown in Fig. [3] 
by the arrow at CEP. From directions not parallel to the tangent line 7 ~ 0.64, as was found in Ref. (23| . 

Below we present in details how the direction corresponding to the Ising temperature axis was found. According 
to the Landau-Ginzburg type analysis of Ref. jl^l, the critical region of the CEP is splited into two scaling regions. 
One of them contains path which are asymptotically parallel to the tangent line of the cross-over curve at the CEP 
along which 7=1 (mean field exponent), while the other one is the complement region of the first where 7 = 2/3. 
In accordance with these, we found that by choosing a path which is near the asymptotically parallel direction, the 
chiral susceptibility as function of the distance from the CEP consists of two straight lines in a log-scale plot having 
slopes 1 and 2/3, respectively. Moreover, as we approach the direction of the tangent fine to the critical curve at the 
CEP, the range of the line with slope 2/3 shrinks and the range of the line with slope 1 grows, as one can see in Fig.[4l 
Asymptotically parallel to the tangent line there is no more crossing of scaling regions, the path remains in the region 
characterized by 7 = 1. 

IV. CONCLUSION 

A resummed perturbation theory was used for the parametrization of the SU{3)l x SU{3)r chiral quark model 
with one-loop accuracy using as input the physical pseudoscalar sector and the masses of the constituent quarks. We 
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found that in a large range of the bosonic and fermionic renormalization scales, within the parametrization scheme 
used the model gives I''* order transition on the T = axis and crossover type transition on the fiB = axis of 
the fJ-B — T-plane at the physical point of the — rriK mass-plane. The renormalization scales were fixed by the 
requirement of minimal deviation of the predicted mass spectrum from the physical one and the critical end point was 
located at /is.CEP = 895.38 MeV and Tcep = 74.83 MeV. The location and the shape of the critical region of the CEP 
were confronted with other results from the literature. We consider that in an effective model a finite answer for the 
actual location of the CEP can be given only by supplementing the chiral nature of the transition with other aspects 
of a true deconfinement transition, like the influence of the conflnement, which can be taken into account through 
the Polyakov loop potential and the dependence of the axial anomaly on the temperature and density. For this more 
detailed knowledge of these phenomena is needed. 

We also studied the surface of the second order phase transition points in the — rriK — fJ-B space using chiral 
perturbation theory for meson and baryons in order to continue the model's parameters from their value at the physical 
point. This surface rises from the boundary of the first order phase transition region on the — rriK plane as the 
baryonic chemical potential is increased. We found that the surface has positive curvature bending monotonically 
towards the physical point of the mass-plane. 
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APPENDIX A: CALCULATION IN THE MIXING SECTOR 



Not all of the scalar and pseudoscalar fields of the Lagrangian are mass eigenstates, there is a mixing between 
the singlet and the octet states. Since physically only mass eigenstates can propagate we have to diagonalize the 
propagators encountered when calculating the self-energies. In the mixing, x — y sector, the diagonal matrix G 
containing the tree-level physical propagators are obtained through an orthogonal transformation performed with 



cos ( 



the matrix O = ^ "sin 6* cos6'J (taii26' — 2m^y/(m^^ — rUyy)) as G ~ OGO^ . Following the convention of 

Ref. (see Eqs. (18) and (20) ) we define G5 — diag(Gcr, G/^) (Gp = diag(G,,', G,,)) in the (pseudo) scalar 
sector. It is easy to show, that one has sa\29p = 2m,-ixy/ J (m'^^.^j. — m'^^yy)'^ + 4:Tnf-i^.y, cos^flp = [1 + (tj^^^ — 



(''^vxx " '^qyyf' + 4TO^^y]/2 in the pseudoscalar sector and sin205 = —Irriaxyl \J {fnlxx ^ ^''-ayy)'^ + ^w^i^yi 

cos^ 6*5 = [1 — (m^ra; ~ "^"ayy) / \J {''^^axx ~ '"^'ayyY + 4'Ti^^y]/2 in the scalar sector. These expressions are to be used 
when expressing the xx, xy and yy components of the internal propagators of Feynman diagrams. In the pseudoscalar 
sector one has 

Gf^ = G^'cos^e'p + G^sin2 6'p, G^^ = G^/ sin^ 6'p + G,, cos^ 6ip, Gf^ = (G^- - G,,) sin ^p cos 6'p, (Al) 

and similarly for the scalar sector with the replacements Gn' G^, G,, — > G/„ and 9p 63. By doing this in the 
non-mixing sector we obtain the self-energies in terms of the physical propagators, while in the mixing sectors we obtain 
the XX, xy and yy components of the self-energies. In this latter sector an additional diagonaHzation is required to 
obtain the self-energies for the mass eigenstates 77', 77, /o and a. 



APPENDIX B: INTEGRALS 



The one-loop fermionic integrals encountered in the calculation are the tadpole Tp{mq) and the bubble Ip{p,mq). 
In the imaginary-time formaHsm the tadpole integral for a quark with mass niq is defined as: 

(Pk 



= -4mq^T J2 / .2^^3 , ^ =■ -4m,T^(m,), (Bl) 



m 

where a;„ = (2n + l)z7rT. After renormalization at finite temperature one has 

^2 ^^2 



1 



^^(^9) = Y^ln-^ + ^ / duj^luj^-mlfpiu:), (B2) 
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where in terms of the Fermi-Dirac distributions for quarks and antiquarks f${i^) = l/(exp(/3(a; =F y^tg)) + 1), fri^) is 

given by fri^^) = —5 [fpi^) + fpi^)] ^iid fiq = //_b/3. (The bosonic tadpole Tglrn) looks exactly the same but with 
friuj) replaced by the Bose-Einstein distribution /_b(w) = l/(exp(/3u;) — 1)). 

There are two types of bubble integrals, the first occurs in the pseudoscalar self-energies and is used in the process 
of the parametrization while the other gives contribution to the scalar self-energies when calculating the one-loop mass 
of the scalars. Both are needed only at zero temperature. Their expressions stand here: 



^p/s{P,'mi,m2) = ^Tr 



d'^k iTp/s{l/: + rni)Vp/s{i>- ^ + m2) 
(27r)4 {k^ ~ ml){{k - pY ~ ml) 



9f [{p^ - (™i T m2)'^)I{p, mi, 1712) - Tpimi) - TF{m2)\ , 



(B3) 



where for pseudoscalars Tp = 75 and for scalars Ts 
a bosonic bubble integral: 



-i. I{p,mi,m2) is the expression encountered when calculating 



/(p,TOi,m2) 



1 



G 



167r2 



2 



arctaii ■ 



, mi 1 , 

■ +ml— p'^ 



■ml 



ml 



G 



2 2 
— mS mf 

-^In^ 

3^ mi 



G 



p^-G 



— inQ{p^ — {mi + 7712)^), for > (mi + 7712)^, < (mi — m2)^ 



(B4) 



arctan ■ 



G 



for (mi — m2)^ < p^ < (m,i + m2)^ 



where G = | (p^ — (mi + m2)^) {p^ — (mi — m2)^) 
easily obtained. 



1 1/2 



For two equal masses the corresponding expression can be 
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